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Abstract 

In this paper an updated Geant-4 simulation toolkit is presented de- 
scribing the interactions in matter of heavy hadrons containing long-lived 
squarks and gluinos. Generic signatures are derived that are applicable at 
the Large Hadron Collider, and the problem of theoretical uncertainties 
arising from model dependence is addressed. For a more detailed descrip- 
tion of the assumptions and techniques used, the reader is referred to ref. 



1 Introduction 

Stable massive particlefS (SMPs) are a feature of many models of physics be- 
yond the Standard Model [H [3] . Exclusively weakly interacting SMPs may be 
searched for utilising the induced missing transverse energy in events where 
they are created. Electromagnetically charged SMPs (e.g. sleptons) may be 
recognized by their ionisation properties as they traverse a detector as well as 
(depending on their mass) their time of flight. Particles carrying colour charge 
on the other hand exhibit a phenomenology of somewhat larger complexity. 
Squarks and gluinos for instance will form hadrons (hereafter termed R-hadrons 
in accordance with litterature) that will then undergo not only electromagnetic 
but also hadronic processes in the traversed matter. This work describes an 
implementation of a Geant-4[3] toolkit describing scattering of R-hadrons in 
matter. The toolkit is described in detail in refs. [T1[S] and is based on previous 
works [3 [7]. 



2 Scattering of heavy hadrons 

When describing the scattering of the heavy hadrons formed by long-lived 
squarks or gluinos as they traverse a detector, two questions become important: 

^Here taken to mean that the decay length is larger than a typical collider physics experi- 
ment. 
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What are the stable hadronic states, and what are their scattering properties? 
Answers to both questions, are of course highly model dependent. Previous 
approaches [3 H] have employed liberal assumptions with respect to the pos- 
sible stable states allowing for all combinations of u and d quarks to combine 
with the heavy sparticle to form hadrons. For this work, a number of different 
mass calculations have been taken into account [HI IB UHl IHl HI] ■ From these 
calculations a consistent set of masses can be chosen as shown in table [1] 



Heavy Parton 


States 


Mass (GeV) 


Squark 


qu, qd 


rriq + 0.3 




qud 


ruq + 0.7 


Gluino 


gqq,gud, gdu, gg 


rUg + 0.7 




guds 


rrig + 0.7 



Table 1: Assumed stable hadrons formed from squarks [q) and gluinos (g), 
together with mass estimates used in this work. The neutral gluino states 
containing a mixture uu and dd pairs are generically denoted by gqq. 

Of particular interest in table [T] is the observation that especially the baryonic 
states are constrained. For instance there can be no tww++ or guuu^^ states. 

Turning to the scattering of the particles themselves, the toolkit used was based 
on that presented in [7]. The software framework and the treatment of nuclear 
effects such as nuclear evaporation and creation of black track particles was 
left unchanged. The innovations were the implementation of a theory driven 
scattering cross section and a different approach to the transmutation of the 
R-hadrons themselves. 

The cross section calculation used stems from the model presented in ref. [13] 
(the Regge model). It is depicted in figure [1] together with the predictions from 
the previous model (refs. [6l[7], hereafter referred to as the generic model). 
For the specific formulae the reader is referred to ref. [1]. Final states in the 
individual hadronic interaction are chosen randomly among the kinematically 
allowed processes with the constraint that meson states should turn into baryon 
states with a 10% probability. Once an i?-hadron becomes a baryon it cannot 
turn back into a meson by subsequent hadronic collisions as was first noted in 

i- 

Another process that must be taken into account is the possibility that neutral 
squark mesoninos may oscillate into their anti-particle. As mixing of such states 
would be highly model dependent we allow only for the cases of no mixing 
and maximal mixing, where there is a 50% probability for any netral squark 
mesonino to turn into its anti-particle. These choices correspond to oscillation 
lengths of zero and infinity, respectively. 
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Figure 1: Cross section per nucleon for the interaction of stop-based and gluino- 
based i?-hadrons. The predictions are shown for the Regge model (sohd hnes) 
and the generic model (dashed lines). 



3 Signatures deriving from hadronic interactions 



The heavy sparticle is typically orders of magnitude heavier than the light 
quark system (LQS) of the R-hadron. As the scale (i.e. geometric cross sec- 
tion) of a wave function is inversely proportional to the square of the mass 
[5], any hadronic interaction is assumed to take place between a nucleon and 
the LQS of the R-hadron. The kinetic energy of the LQS system scaling with 
Mlqs Ekin.Totai therefore represents a measure of the available energy in any 



collision. This is consistent with the observed energy loss per hadronic inter- 
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Figure 2: Left: Energy loss per hadronic interaction for stop hadrons in iron 
is shown for the Regge model with and without mixing as well as the generic 
model. Right: Total energy loss for R-hadrons traversing 2 m of iron 

action shown in figure [2] (left). Generating gluino pair-production events with 
Pythia[T3j for pp collisions at a centre-of-mass energy of 14 TeV, kinematic dis- 
tributions were obtained for R-hadrons. R-hadrons were subsequentially simu- 
lated traversing 2 m of iron to estimate their behaviour in an LHC experiment 
[15] . The results showed that the different types of R-hadrons would undergo 
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~3-8 hadronic interactions and deposit energy as shown in figure [5] (right). 
These observations imply that R-hadrons at the LHC experiments will inter- 
act hadronically but that they will not through these interactions give rise to 
anything that might be interpreted as hard jets. 

Using similarly generated kinematics as for the energy loss observations, the 
distributions of R-hadrons on the different flavours in table [1] were studied as 
a function of penetration depth in iron. As is shown in figure [3] (left), gluino 




Figure 3: Left: Flavour distribution of 300 GeV gluino R-hadrons as a func- 
tion of penetration depth in iron. Right: Percentages of different charges and 
stopping fractions at a penetration depth of 2 m in iron for 300 GeV R-hadrons 

R-hadrons have a ~50% probability of being baryons after traversing 2 m of 
iron. Of the remaining R-mcsons, roughly half will be electrically neutral. The 
stopped fraction of the R-hadrons is on the order of a few percent. Comparing 
these numbers to the stop and sbottom cases, one needs to consider both squarks 
and anti-squarks for the minimal and maximal mixing scenarios. This is done 
in the table in figure [3] (right) where the gluino numbers are also included for 
comparison. As can be seen, the fraction of stopped R-hadrons is in all cases 3- 
8%. Gluino and sbottom hadrons in all scenarios are found to be predominantly 
neutral while the stop is found to have a charged fraction of roughly 2/3. These 
observations are not greatly diluted when mixing is enabled. 

It turns out that there is little impact on these conclusions by varying the 
meson/baryon mass splittings in the model. Varying the mass splitting by 1 
GeV moves the predictions less than one percent. Varying the hadronic cross 
section by a factor of 2 produces an effect of 0(10%) in the predictions. The 
scattering cross section thus remains a key ingredient in the estimation of the 
theoretical uncertainties for the scattering of heavy hadrons containing coloured 
sparticles. Should it turn out that more baryonic states of squark or gluino 
hadrons exists, the prediction of this model would of course also be altered. 
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